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Thermodynamics of the spin | XXZ model is studied in the critical regime using the 
quantum transfer matrix (QTM) approach. We find functional relations indexed by the 
Takahashi-Suzuki numbers among the fusion hierarchy of the QTM's (T-system) and their 
certain combinations (F-system). By investigating analyticity of the latter, we derive a 
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closed set of non-linear integral equations which characterize the free energy and the corre- 
lation lengths for both (c^cr~) and (o~jCrf) at any finite temperatures. Concerning the free 
energy, they exactly coincide with Takahashi-Suzuki's TBA equations based on the string 
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• hypothesis. By solving the integral equations numerically the correlation lengths are deter 

, mined, which agrees with the earlier results in the low temperature limit. 
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1 Introduction 

X. 

In this paper we study the spin | XXZ model at finite temperature based on the recently 



developed quantum transfer matrix (QTM) approach |§-|17]. We shall deal with the "root of 



unity" case in the gap- less regime. Namely, the anisotropy parameter has the form A = cos^ 



with pq any rational number not less than 2. (See fl2.ip .) We derive the non- linear integral 
equations that characterize the free energy and the correlation lengths for both ((rfcf) and 
(cjcr?) at any finite temperatures. 

Thermodynamics of the XXZ model is a classical and by no means fresh problem at least 
as far as the free energy is concerned. It goes back to 1972 that Takahashi and Suzuki [p~8[ 
took the thermodynamic Bethe ansatz (TBA) approach [19] to the free energy based on the 



elaborate string hypothesis. They selected, as allowed lengths of strings, a special sequence of 
integers nj which we call the Takahashi-Suzuki (TS) numbers. The resulting free energy yields 
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correct physical behaviours in many respects. Actually this is one of the best known example 
among many successful applications of the TBA and string hypotheses. However there is also 



some controversy in string hypotheses themselves [gy, 21, ^2|, in view of which those successes 
are rather mysterious. 

This is one of our motivations to revisit the XXZ model with the recent QTM method. It 
integrates many ideas in the statistical mechanics and solvable models |l|]-[]l7| and has a number 
of advantages over the traditional TBA approach. It only relies on certain analyticity of the 
QTM, which can easily be confirmed much more convincingly by numerics. Moreover it enables 
us to systematically calculate the correlation lengths beyond the free energy for a wide range 
of temperatures. See [12] for A > 1 case. Roughly, the QTM method goes as follows. First 
one transforms the ID quantum system into an integrable 2D classical system based on the 
general equivalence theorem jj], ^]. The QTM T\ is a transfer matrix propagating in the cross 
channel of the latter. Despite that the original ID Hamiltonian H is critical, the QTM T\ can 
be made to have a gap. Therefore the formidable sum Tre _/3 ^ (/3 = 1/kgT : T is temperature) 
can be expressed as its single eigenvalue which is largest in the magnitude. Furthermore the 
correlation lengths are obtained from the ratio of the largest and the sub-leading eigenvalues 
of T\. To evaluate them actually, one must however recognize a price to pay; now the QTM 
T\{u) itself becomes dependent on the fictitious Trotter number N through its size and also a 
coupling constant as T±(u = const ^) [^-fllPf]. This makes it difficult to determine the spectra 
of T\ by a naive numerical extrapolation as N —> oo. A crucial idea to overcome this is to equip 
the QTM with another variable v and to exploit the Yang-Baxter integrability with respect to 
it; [Ti(u,v),Ti(u,v')] = [11]. Here u and v play the role of the (inverse) temperature and 
the spectral parameter, respectively. Furthermore one introduces some auxiliary functions of 
v, which should realize somewhat miraculous features. Their appropriate combinations should 
have a nice analyticity that encodes the information on the Bethe ansatz roots of Ti(u, v). Once 
this is achieved, one can derive a non-linear integral equation which efficiently determines the 
sought eigenvalues of T\. The Trotter limit N —* oo can thereby be taken analytically. The 
most essential step in this method is to invent such auxiliary functions and their appropriate 
combinations. There are some interesting variety of choices for them in various models |l2]| — p^] . 

Back to the XXZ model our finding is that such auxiliary functions can be given by the 
QTM's {T n .-i\rij : TS numbers}, which is the subset of the known fusion hierarchy of com- 
muting transfer matrices whose dimensions of the auxiliary spaces are precisely the TS numbers 
nj. We will show that {T nj ._i} satisfy functional relations among themselves (T-system) and 
so do their certain ratios {Yj} an elaborate one (Y-system). See (|4.6|)-( fi~9| ). Especially there is 
a special identity (fO^) among {T nj _i} that holds only at rational values of po and makes the 
Y-system close finitely. Besides the peculiarity at general roots of unity, such use of the fusion 
hierarchy as the auxiliary functions originates in the studies of finite size corrections [23, p4| [. 

As for the free energy we thus obtain the integral equations identical with Takahashi-Suzuki's 
TBA equations but totally independently of their string hypothesis. We shall further study the 
second and the third largest eigenvalues of T\ for po integer. They are related to the correlation 
length £ of (crTcrr) and (a^af), respectively. In contrast with the largest eigenvalue, now the 
zeros of the T nj —\ come into the "physical strip" spoiling the nice analyticity. Nevertheless 
we manage to identify their patterns and derive the "excited state TBA equations". Solving 
them numerically we determine the curve £ = £(/?). Especially the low temperature asymp- 
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totics lim^oo £((3)/ (3 agrees with the known result |T^, ^] with high accuracy for the both 
correlations. 

Our formulation here using the 2 variable QTM Ti(u, v) and fusion hierarchies is based on 
[ p| , 16]. There are similar approaches in the context of integrable QFT's in a finite volume 



27, 28|. 



It has been known for some time that solutions of Y-systems can curiously be constructed 
from T-systems [23, 29]. By now this connection has been generalized to arbitrary non- twisted 
affine algebra Xr ioi the associated Y-system |3(| and the T-system [p9|] . (See also fl3lf .) In 



this sense, our results here display a further connection of such sort for U q {sl{2)) at q general 
root of unity. 

The layout of the paper is as follows. In section 2 we formulate the XXZ model at finite 
temperature in terms of the QTM T\, In section 3 we give the fusion hierarchy {T n _i} of 
QTM's and their eigenvalues. A functional relation (T-system) valid for general po is also 
given. In section 4 we construct the Y-system out of the T-system. The former closes finitely 
due to the special functional relation (|4.2| ) valid only for rational p$. In section 5 we derive 
integral equations for the free energy, and in section 6 for the correlation lengths of (a^a^) 
and (crjof). Section 7 is a discussion. Appendix A recalls the definition of the TS numbers and 
the related data. Appendices B and C contain a check of the analyticity of the Y-functions for 
the free energy. Appendix D is devoted to the free fermion case po = 2, which needs a separate 
treatment. 



2 Quantum transfer matrix 

The Hamiltonian of spin | one dimensional XXZ model on a periodic lattice with L sites is 



H 



H 



33+1 



3=1 
J 
1 



(2.1) 



Here o"J,(jJ,cr| are the local spin operators (Pauli matrices) at the j-th lattice site and J 
is a real coupling constant. We shall consider the model with the anisotropy parameter A 
in the critical region — 1 < A < 1. Due to the invariance of the spectrum of (2.1) under 
the transformation (J, A) — > (— J, — A), we can further restrict the range to < A < 1 and 
introduce the parametrization: 



A 

Po 



cos6 0<9< 
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2' 
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Po > 2. 



(2.2) 
(2.3) 



ni/po 



The model is associated with the quantum group U q (sl(2)) at q = e 

In order to consider its thermodynamics we relate it to the six vertex model. This is a two 
dimensional classical system whose Boltzmann weights are given by 
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where 



sin 2^ 



(2.5) 



Let V be a two dimensional irreducible module over U q {sl{2)). As is well known the quantum 
ii-matrix R € End(V ® V) with the above matrix elements and the spectral parameter v 
satisfies the Yang-Baxter equation (YBE) (cf.( |3.3| )). To relate the six vertex model with the 
XXZ model, consider a two dimensional square lattice with N rows and L columns. We shall 
assume that ./V is even throughout. Define the (auxiliary) transfer matrix T A (u,v) as 



T A (u,v) = Tv 



Vi®V 2 
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(2.6) 



See also Fig.[l]. Here and in what follows V\ 



V N = V{ 



V[ = V. Operators 



diagrammatically shown as in (2.6) are always assumed to act on the states in the bottom line 
to transfer them into those in the upper line. Using the identity (Pjj+i is the permutation 
operator acting on Vj <g> Vj+i), 



v=0 



P 3 3+lj- R 3 3+l( V ) 



J sin 



H 



v=0 

i.i+i + 2 A 



we expand Ta(u,0) as 



T A (u,0) = 1 + 



26u 



H + 



JLA\ 



\+0(u z 



(2.7) 



Jsin6» V" 2 

This formula represents as an equivalence of the XXZ model and the six vertex model. In fact 
we can go further to the finite temperature case. From ( |2.7| ) it follows that 

(3 J sin 9 



-P[H + ^-a)) = lim T a (u n ,0)t. 

A II N— >oo 



exp ( • i // — ^ ; ; i , / i ! " \ ■ u i - ■ «.\ 
Thus the free energy per site / of the XXZ model is given by 



ON 



f = - lim lim ^— In (Ty v ® l T a {u n , 0) a ) - ^A. 



(2. 



(2.9) 



However, eigenvalues of the transfer matrix T a (un,Q) are infinitely degenerate in the limit 
mat ^^5° 0, therefore the the trace renders a serious problem. To avoid this, we introduce the 
following trick; we rotate the lattice by 90° (see Figjl|) and rewrite the free energy as 



/ = - lim lim — In ( Tr y ®;vTi(iiAr, 0) J 

L^ooN^ooLp \ 



—A, 
2 ' 



(2.10) 



where 



Ti (u, v) = Try/ 
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Figure 1: Relation between Ta(u, v) and Ti(u,v). 



We call Ti(u,v) the quantum transfer matrix (QTM). Note that due to the YBE ( |3,3[ ), the 
QTM is commutative as long as the u variable is taken same: 



[T 1 (u,v),T 1 (u,v')]=0. 



From now on, we write the /c-th largest eigenvalue of the matrix T\{u,v) as T^ k \u,v). Since 
the two limits are exchangeable as proved in B H, we take the limit L — > oo first. Noting that 



there is a finite gap between lim^v- 

/ = 



(1) 



(2) 

(utv,0) and limTv^oo^i (i*7V,0), we have 



J 



- lim lnT«(^,0)-^A. 



(2.12) 



Namely, the problem of describing the thermodynamics of one dimensional quantum systems 
reduces to finding the largest eigenvalue of the QTM of two dimensional finite systems (to be 
exact, finite in the vertical direction only). Equation ( 2.12| ) is a finite temperature extension of 
the equation ( |2.7| ). 

In this approach, the thermodynamical completeness — lim^o (5f = In 2 follows easily from 
(0,0) = 2, which is obvious from i?i2(0) = P\2- (As for the combinatorial completeness see 
fffil including the higher spin cases.) 

Most significantly this method makes it possible to calculate some correlation length (k > 
2) at finite temperature. To see this let Qi = 1<2>- • -(8)Q(8)- • -<8>1 6 End(V r(glL ) be a local operator 
acting on the i-th site V( via Q = J2ii'=±i Qe',e^e',e- Here Egr^ denotes the 2 by 2 elementary 
matrix and Qg g is the matrix element. (a z = En — -E_i 5 _i,o" + = Ei—%,a~ = E-n.) Given 
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Q we introduce the operator S(Q\u, v) <G End(V® N ) by 



S(Q\u, v) 



E Q 

ei'=±i 
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Then for the local operators Qi,Pj E ~End(V® L ) their finite temperature correlation function is 
expressed as 

Ti v ®l PjQi exp(-(3H) 



(PjQn 



lim 



L^oo Tiy<s,L exp(—f3H) 

Tr v ^ N (s{P\u N ,0)T 1 (u N ,Oy- i - 1 S{Q\u N ,0)T 1 {u N ,0) L+i ^- 1 ) 
= lim lim — 

JV—ooL^oo Tt v ®nTi(un, 0) L 

where we have exchanged the two limits. Suppose that and are the operators 

such that the matrix elements of S(Q^ \un, 0) and S(Q^\un, 0) between the eigenspaces 
for T^\un, 0) and X^ fe ^(^,0) are zero for 1 < kl < k and nonzero for kl = k. Then setting 
(P, Q) = (Q( fc )t 5 QW) i n the above, we have 

j-i 



(Q^Q\ k) ) 



lim 



Tl k \ UN ,oy 



Fitting this with exp(— (j — i)/£k) hi the limit j » i we have 



k > 2, j > i. 



— = — lim In 



Ti 1] (u N ,0) 



(2.13) 



As seen in section 6, £2 and £3 are the correlation lengths of (<7j~<7j ) (Q^ = a ) and (crjof) 
(Q( 3 ) = a z ), respectively. 

3 T- system 

To study Ti(u, v), we embed it into a more general family of transfer matrices and explore the 
functional relations that govern the total system. For this we consider the fusion hierarchy 
{T n _i(u, v)} defined by 



T n - 1 (u,v)=Tt\ 
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where the trace is over the n dimensional irreducible auxiliary space depicted by the dotted 
line. To be explicit, we give the constituent fusion Boltzmann weights. 
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(3.2) 



Here j,f £ {1, . . . ,n} and \j - j'\ = 1. 

e ii'' e jj' are arbitrary parameters such that Cjj'Cj'j — 
e'jjiCjij = 1- If they are 1, the six vertex case n = 2 of these weights reduce to (^4|) under 
the identification of j = 1 and 2 states with +1 and —1, respectively. The Boltzmann weights 
satisfy the YBE: 






(3.3) 



From the picture ( |3.l|) one sees that the members of the fusion hierarchy are all commutative 
for the same u: 

[T n _i(u,v),T„>_i(u,?/)] = 0, 

due to the i?-matrix intertwining the n and n' dimensional representations. Thus they can be 
simultaneously diagonalized and the eigenvalues (also written as T n ^\(u, v)) are readily obtained 
in the dressed vacuum form: 



T n _i(it, v) = <t>(v - i{u + n + 2 - 2j))<f>(v + i{u-n + 2j)) 
J'=l 

Q(v + in)Q(v — in) 
X Q(v + i(2j - n))Q(v + i(2j - n - 2)) ' 



(3.4) 



sm i 

m 



Q(v) =nsh-(t;-w i ). (3.6) 

i=i 

Here m £ {0, 1, ... , -/V/2} is the quantum number counting the (— l)-states on odd sites and 
(+l)-states on even sites. The dressed vacuum form is built upon the pseudo vacuum state 
((+1) <8> (—1)) 2 , which corresponds to m = 0. {tOj} is a solution of the Bethe ansatz equation 
(BAE): 

N_ 

I shf(u;j + i(u + 2))shf (uj - iu)\ 2 _ Q(ujj + 2i) 



shf^- - i(u + 2))sh|(w J - + m) / - 2i) ' 



(3.7) 
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The largest eigenvalue of T\(u,v) lies in the sector m = N/2. Note that T_i(it, v) = and 
Tq(u, v) = <p(v — i(u + l))<p(v + i(u + 1)). An important property is the periodicity 



T n -i(u,v) = T n -i(u,v + 2p i). (3.8) 

Let us present the functional relations among the fusion hierarchy. For any v € C and 
integers n > y > 1, the following is valid, which we call the T-system. 

T n -i(v + iy)T n _i(u - iy) = T„ +y _i(?;)T n ^_i(?;) + T y -i(v + in)Ty-\{v - in). (3.9) 

Hereafter we shall often omit the common u variable to simplify the notation. The proof of 



this equation is direct by using the expression (3^4). Representation theoretically, it is a simple 



consequence of the general exact sequence in [33[ as explained in p9] for y = 1. In general the 



T-system ( [3.9D extends over infinitely many transfer matrices. However, as we shall see in the 
next section, for rational po there is a special functional relation (|4.2[ ) that makes the associated 
y-system closes finitely. 

4 y-system at Root of Unity 

From now on, we shall concentrate on the case when p$ > 2 is a rational number and treat the 
free fermion case po=2 separately in Appendix D. Consider the continued fraction expansion of 
Po 

Po = ^i + 1 , (4.1) 

V2 H 



1 



which specifies a > 1 and ui,...,v a £ Z>i. From the assumption po > 2, we have u± > 2. In 
fact V\ = 2 is allowed only if a > 2, and po = v\ > 3 is assumed if a = 1. 

In Appendix A we recall the sequences of numbers {mj}^ 1 , fe}?=o > {yj}j=-ii { z j}j=-i 
and {%}j>i introduced in |18|| . The last one is the TS numbers. We shall also introduce its 
slight rearrangement {«j}j>i and a similar sequence {wj}j>i related to the "parity" of the 
TS strings. They are all specified uniquely from pq. With those definitions we now describe a 
functional relation of the T n _i(u), which is valid only at the root of unity and is relevant to our 
subsequent argument. 

Ty a +y a _ 1 -i(v) = Ty^y^iv) + 2(-l) mz <*T ya _ 1 _ 1 (v + iy a ). (4.2) 

Here m is the number of the BAE roots in (|3.6| ). The proof is straightforward by using (3.4) 



Lq ) and Q(v + 2iy a ) = (— l) mZa Q(v). When a = 1 hence po = u i, ( f4.2|) reduces to a simple 
relation r^(u) = T vi _ 2 (v) + 2(-l) m T (w + ii^). 

For 1 < j < jmax '■= ro 8 -l, let < r < a— 1 be the unique integer satisfying m r < j < m r+ i. 

Set 

y/x + Wjp )T fij . +1 -y r -i(z; + Wjpp) 

j r yr _i(u + inj+i + ivjjPo)T Vr -i(v - in j+ i + iwjp ) ' 

rn j+1 -i(v + iy r + »WjPo)T fij . +1 _i(T; - zy r + iwjpo) 
Ii, r _i(u + zn i+ i + wojPo)-T yr _i(t; - + tWjp ) 
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K(V) = " 3ma * - JrnavWJ ^ 

T y a -i~l{V + Wa + lWj max PO) 

where in (|45|) rij max = y a — Va-x and Wj max = z a — z a -\ — 1 in accordance with flA.lCf) and 
(|A.11P . We also set Y Q (v) = and Y-i(v) = oo. Thanks to the T-system and 



( [4.4|) are equivalent. We find that {Yj(t>)}j=i 1 and K(v) close among the following finite set 
of functional relations, which we call the Y-system. 

Theorem 1 

For m r -i < j < m r — 2 (1 < r < a), 

Y 3 {v + ip^Yjiv - ip r ) = (1 + Yj^v)) 1 - 25 ^-! (i + Y j+1 (v)), (4.6) 
for j = m r — 1 (1 < r < a — 1), 

Yj(v + ip r + ip r+ i)Yj(v + ip r - ip r+1 )Yj(v - ip r + ip r+1 )Yj(v - ip r - ip r+ i) 

1 — 25i„ r 

(1 + Yj-i(v + ip r +l))(l + *J-l(« - iPr+l)) | (1 + + *Pr)) 

x(l + Y i+1 (« - ip r ))(l + + - ip r +l))(l + Yj(v ~ + iPr+l)), (4.7) 

i + y ma _ 1 W = (i + ^)) 2 , (4.8) 

K(v + - %) = 1 + Y ma „ 2 (v). (4.9) 



This can be proved by combining the T-systems (|3.9| ), ( [4.2| ) with the definitions of 
{toj}, {pj}, {yj}, {nj} and {^j} in Appendix A. When u r = 1, ( ^4.6|) is void and ( f4.7|) holds 



for j = m r — 1 = m r _i. See (A.l). 



In the case po £ Z>3, the Y-system has a simple form (1 < j < po — 2) 

Yj(v + i)Yj(v - i) = (1 + Yj-iMX 1 + ( 4 - 10 ) 

l + y po _ 1 ( U ) = (l + ^(t;)) 2 , (4.11) 

+ - i) = i + r P0 _ 2 («), (4.12) 



where 



T (v + i(j + l))T (v-i(j + 1))' 
Tj(v + i)Tj(v - i) 
T (v + i(j + l))T (v-i(j + 1)) ' 



1 + Yj (v) = ^^ 7 ^f i V^J).. 1 I i~v\" ; ( 4 - 14 ) 



= (-r r ;-; M (4.i5) 

for 1 < j < Po — 2. Due to the property of the TS-number, Yi(v) and 1 + Y\(v) are always given 
by setting j = 1 in (4.13) and ( |4.14j ) for arbitrary pq > 2. 



A similar Y-system has been considered in [|3j] in a perturbed CFT context. 

5 Integral equation for free energy 

Let \k) be the eigenvector corresponding to the k-th largest eigenvalue Ti k \u,v) of T\(u, v). 
We define the k-th (not necessarily fc-th largest) eigenvalue T^\(u, v) of the auxiliary QTM 
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T n -i(u,v) by r n _i(«,u)|Jfe) = T^\(u,v)\k). Let {Yf>} (and be the Y-functions con- 

structed from as in ( fO| ) - ( |4.5| ). In this section, we study the analyticity of {Yj (u, v)} 

and K^(u, v) in the complex u-plane. Then we derive the integral equations which characterize 
the free energy. 

An advantage in the present approach lies in the fact that the analytic assumption given 
below can be explicitly checked numerically keeping the Trotter number N finite. We have 
performed numerical studies with various values of po, (3 and N in determining the location 
zeros of fusion QTM's. For example, the zeros for T^]_i(u, v) for po = u = ±0.01, n=2,3,4,5 
N = 16 and N = 32 are plotted in Fig [2L Guided by them we have the following for u small 
(typically \v\ ~ 0.01). 

Conjecture 1 All the zeros of T^\(u,v) are located on the almost straight lines = 
±ra mod 2po- 

This coincides with the observation in the XXX model if one forgets the periodicity 2po in the 
imaginary direction. It is also consistent with the functional relation ( |3.9| ). The deviation from 
the straight line is very small (~ 10~ 2 at most) as seen in the figures. It becomes smaller as 
u —* 0. Once Conjecture 1 is assumed, we can identify the strips in the complex f-plane in 
which Yj(v) or 1 + Yj (y) are analytic, nonzero and have constant asymptotics at v = ±oo. 

We call this property ANZC. In Appendix B we verify that Yj^fv), for example, is ANZC in the 
strip < x whenever the combination Yj(v + ix)Yj 1 \v — ix) takes place in the Y-system 
(|4,6D-( |4~9| ). Apart from the exceptional Case 1, 2 and 3 listed below, this makes it possible to 
transform most of the Y-system into integral equations defined on the real axis quite easily. 
This is a consequence of a simple lemma. To present it we let S[x] denote the strip E [— x, x] 
in the complex v-plane {x € R->o)- Then we have 

Lemma 1 Suppose the functions gi(v) satisfy 

g (v - iv )go(v + iv ) = ]J 9j(v - ivj)gj{v + ivj), (5.1) 

where Vj > are real numbers and vq > Vj(j > 1). Assume further that gj(v) is ANZC in the 
strip S[wj] for some Wj > Vj for j > 0. Then the above functional relation can be transformed 
into the integral equation 

/oo 
Rj(v — v') In gj(v')dv' + constant, 

RAv) = — r e lkv ^4dk, 

where the constant is determined by the asymptotic values of the both sides. 

The proof uses Cauchy's theorem and the fact that the ANZC function gj(v) admits the 
Fourier transformation of its logarithmic derivative. 

There are few exceptions to which the above lemma can not be applied directly: 

• Case 1. j = 1 in Q, 

• Case 2. r = 1 when v\ = 3, u% = 1 and V3 ^> 1 in ( [4.7[ ), 
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Case 3. r = 1 when v\ = 2 in (|4~ 



Nevertheless, they can still be converted into integral equations after a suitable recipe. Let us 
explain this for the most important Case 1 below. 
Case 1 in (4.6) is explicitly given by 

Y} 1] (v + i)Y} 1] {v -i) = l + Y^\v). (5.2) 

Yi(v) possesses zeros (resp. poles) of order N/2 at ±(1 + u)i (resp. ±(1 — u)i) mod 2p$i in 
the strip <S[1] for J > (resp. J < 0). (Note that u = ujy is a small quantity given in Q2.8| .) 
Thus the lhs of ( |5.2p does not meet the condition for Lemma 1. A simple trick, however, makes 
it applicable. Define a modified function 

y(i) ( u v ) = Y^\u,v) 

1 1 ' ^ (thf(u-i(l±u))thf(u + i(l±^)) ±JV / 2 ' { } 

where the + and — signs in front of u and N are chosen according as J > and J < 0, 
respectively. Then Y^ (v) has the ANZC property in «S[1]. Due to the trivial identity th^(x + 
i)th.j(x — i) = 1, Y-y (v rfc i) in the lhs of ( |5.2| ) can be replaced by (v =L i). Now the lemma 
applies. The asymptotic values of both sides can be immediately evaluated from the explicit 
results on the T-functions. Then we have, 

m N 7T 7T 

]nY^'(u,v) = ±— ln(th-(u-t(l±«))th-(t; + i(l±u))) 



+ 



A 1 7r(l> — l/) 

oo 4ch-^ i 



L—inCi + yW)^')^'. (5.4) 



2 



Cases 2 and 3 are discussed in Appendix C. In this way all the Y-system can be transformed 
into coupled integral equations. For finite N one can evaluate Y^'s given by (4.3)-(fO|) and 
( |3.4[ )-( [3^ ) directly from the BAE roots. Or one can solve the integral equations numerically. 
We have checked that the two independent calculations lead to the same result up to ./V = 40. 

Let us proceed to the Trotter limit A — > oo. From now on we write the Y-functions in the 
limit as 

r]f\v) = lim Y. (fc) (ujv,tO, (5.5) 
k<®(v) = lim K {k \u N ,v). (5.6) 

N—*oo 



Apart from the Y-functions the A-dependence enters (5.4) only through the "driving term". 
Its large A limit can be taken analytically as 

N ( 7T 7r \ f3nJ sm.0 



1\ I 7T 7T 

lim — In th— (v — ill ± uu))tb.— (v + ± un)) . 
N^oo 2 V 4 V v 77 4 V v "J 29ch(7rv/2 



We thus arrive at the integral equations for rj^ and which are independent of the fictitious 
Trotter number A. 



ln ^ )(u) = " 2 ln|w2) ^ + (1 " 25 ^-ui)sr * ln(l + ^)(v) 

+ s r * ln(l + Tjj^.i){v) for m r _i < j < m r — 2, j > 1, 1 < r < a, (5-7) 
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+ s r +i *ln(l + for J = m r — 1, 1 < r < a, (5.8) 

In = Sa *ln(l + 7 ? «_ 2 )( U ), (5.9) 

where j4 * -B(u) denotes the convolution Jf^ A(v — v')B(v')dv' , and 

1 j f°° Jkv ch(p r -p r+1 )k 



ip r ch^ J-oo 47rch(p r .fc)ch(p r+ ifc) 

The set of the equations closes by one further algebraic equation: f7^_i(«) = k( 1 \v) 2 +2k( 1 \v). 
Under the identification rjj (v) = f]j{v) and k^[v) = k(v), the eqs. ( |5.7[ )-([5T9|) are nothing 
but the TBA equation (3.17) in [jl^] with zero external fieldj^] 

To obtain the free energy per site recall that T±(u, v) satisfies the inversion identity 

T^iv + i):Zi (1) (u - i) = T (v + 2i)T (v - 2i)(l + ^iV))- 

Again, the ANZC property of the both sides leads to 

]nTP(u,v) = si *ln(l + y} 1] )(v) +ln<p(v + i(u + 2))<p(v -*(u + 2)) 

+ N r d 4j 



- dfc_ ito sh*ush(l-f) 



2k" chfeh^ 



Calculating the limit in (2.12) we obtain 



9 

where 



/ = - 27rJ sm - / ai(v)s 1 (v)dv - k B T [ s^v) ln(l + riP(y))dv, 



I sin( 

ax[v) 



2po ch(9v) — cos 9 
This coincides with eq.(3.12) in [lj| under the convention /cb = 1. 

6 Correlation length 



Let us study the correlation lengths of (crT a i ) and (c|of) along the scheme ( p. 13 ). They are 



relevant to the second and the third largest eigenvalues 7\ {u,v) and 7\ (u,v) of the QTM, 
respectively. The former lies in the sector m = N/2 — 1 and the latter inm = N/2, where m 
is the number of the Bethe ansatz roots in ( |3.6D , In this section we shall exclusively consider 
the case po G Z>3 and J > 0, when the Y-system and Y-functions take the simple forms 

(BHED- 

First we need to allocate the zeros of T^Zi{u,v) (k = 2,3) for 2 < n < po in the complex 
v-plane when u is negative small. Based on numerical studies, we have the following for u 
negative small (typically u ~ —0.1). 

1 In their second equation, the range 1 < i < a should be corrected as 1 < i < a. Also in their third equation 
di should be replaced with di. 
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(2) (2) 12) 

Conjecture 2 For 2 < n < po, T„_ 1 (ti, v) has two real zeros ±Q-i f or some Cn-i e R>o- 
AZZ the other zeros o/T^_ 1 (u, w) (2 < n < po) are located on the almost straight lines = 
±n mod 2pQ. 

For example see Fig. |3] showing the zeros of T^^u, t>) for the casepo = 5, u = —0.1, n = 2,3,4,5 
and iV = 20. The main difference from the largest eigenvalue case is the presence of the two real 
zeros for n < pq. Their absence for n = po can be explained as follows. T^\{u, v) in (fO|)-(|3lf) 
is a Laurent polynomial of e^ v . When m = N/2 — 1, its highest/lowest terms are proportional 
to S1 ^jf e^~ v . This is vanishing when n = po, therefore the number of zeros decreases from N 
to N - 2. As a result Y^ ] _ 2 (v) tends to zero as e ' v ' for v — > ±oo. 

As for the third largest eigenvalue we have the following for u negative small (typically 
u ~ -0.1). 

Conjecture 3 T^} l (u, v) has two real zeros ±Ci-i f or some Ci-i £ R>o f or n < Po an d a 

(3) (3) 

double zero at Q -i = f or n = Po- All the other zeros of X) l _ 1 (u,v) (2 < n < po) are located 
on the almost straight lines = ±n mod 2pQ. 

See Fig.|| showing the zeros of T^ i } 1 (u, v) under the same conditions with T^} 1 (u, v). Again the 
main difference from the largest eigenvalue is the two additional zeros on the real axis. 

When v -> ±oo, the Y-functions Yj k) and built from via (Hll) - (p^) have the 
asymptotic values 

/ q n sin 7r ( J+2 ) sin ^ 

(2) ^ P0 w Y p^j(j + 2), (6.1) 

J sm — J 

po 

K {2) -» -1, K {3) -^Pq-1. (6.2) 

To apply Lemma 1 to the Y-system ( 4.10| )-( [4.12| ), we modify the y-functions as 



~ Y (v) 

V^) = - J m7i forl ^^^o-2, (6.3) 



K {k \v) = , (6.4) 

thf( v + c^ 2 )thf( v -c^ 2 ) 



where 



F[ k \v) = {thj(v + i(l + u))thj(v-i{l + u))}^(thjv) 25 ^g {k) (v) forp = 3, 
F[ k \v) = thj(v + 4 k) )th^(v-4 k) ){thj(v + i(l + u))thj(v-i(l + u))}i forp /3, 
Fj k \v) = th j(« + C)5i)tli J(« - Cj5? a )th^(t; + C^th^Ct; - cji) for 2 < j < p - 3, 



4 V 'j+i' 4^ 4 
- -i|(u + c£L 3 )th-(u-C^ 3 

The factor </ fc ) (u) defined by 



^io-2W = th^^ + C^s^h^^-d^)^^) 2 ^^)^) forp ^3. 



9 (k) {v) 



exp (-^thf v) for k = 2 
1 for k = 3 
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has been included to compensate the singularity caused by Y po _2(v) tending to zero as e p o 



I "I 



at v = ±00. 2 The zeros v = of T^'{un, v) depend on N and converge to some finite 
values in the Trotter limit N — > 00. By abuse of notation we shall also write their limit as 
■ (Cpo-l = i s valid irrespective of N.) Proceeding as in the free energy case, we get the 
non-linear integral equations obeyed by rjf^ and kP^: 



In 7]i(v) 

In r][ k \v) 
him (v ) 



where 



' TTV ' 

2 / 



PttJ sin 9 
' 20ch( 

•2^,3 1: 
PttJ sin 9 
' 29ch( 



+ si * In ((1 + K^fh^) (v) + 4,2 (ni - |«th™) 



7T 



+24,3 lnth— v for p = 3, 



+ fll * ln(l + + In \ th^-(v + CDthT^ " (2 



TTV ' 

2 > 



7T 



7T 



Sa * ln(l + r/fj x )(l + v $%)(v) + In jth^ + Cj?i)th^ - cjft)} 
+ i n |thj(v + Cf } i)th^(t; - Cf i)} for 2 < j < p - 3, 
Sl * In ((1 + ^)_ 8 )(1 + kW) 2 ^)) („) + In [th£(* + <Mjfo\{v - Q, 



(6.5) 
for po / 3, 
(6.6) 



(k) > 
Po-3/ 



7T 



+24,3 lnth— i) + <5 fe ,2 7" 



7T , 7Tt> 



-vth 



Po 



for po 7^ 3, 



Sl * ln(l + r/W.aXu) + In |thj(t; + C^ 2 )thJ(« - c£L 2 )} + 4,2™, 



(6.7) 

(6.8) 
(6.9) 



exp ( f ( irth^ 



7 T1J 1 

2 ch4? 



for fc = 2 

^ 1 for A; = 3 

Here the integration constants have been fixed from the asymptotic values (|6.1| ) and ( |6.2| ). 

(k) / ffcJT 

addition to these we need to impose the consistency condition coming from T- (Q - 



In 

0, 



which determines the real zeros {±ej | Q k ' > 0,k = 2,3, j G {1, • • • ,po — 2}}. (Cpo-1 = ^0 
From (gig ) and (§1| this leads to setting T/j fc) (cj fc) ± i) = -1 in (|0|)-(10|). Explicitly they 
read 



(3) 



for po = 3, 



PttJ sin 9 



, _ + Sl * l n ( (1 + K M)* h W) (C W + i} _ f *( C <*> + i)th? 



2#sh( 



+4,3 < 7TZ + In 



sh(fCj 3) )+ i 
sh(fc{ 3) )- 



0. 



for po + 3, 
. flirJ sin 



2<9sh( 



for 2 < j < po - 3, 



^ + s 1 *ln(l + 7?? ) )(cf ) +i)+ln 



sh(fcf fc) )+^h(fcf; 



(6.10) 

+ 7TZ = 0, 

(6.11) 



2 This is a distinct feature of the present case compared with [^, ^7], 
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+ ln — TTk — . : - ^ : + H — iw, — In. \+m = o» (6-12) 



for p / 3, 



•1 * ^ ((1 + C 3 )(l + « (fc) ) 2 fc ( *>) (C 2 + + In 



sh(fcff-2) + fch(Kff- 3 ) 
.sMfC^ 2 )-ich(fcS- 3 ) 



(6.13) 



where the convolutions should be interpreted as 



Here p. v. means the principal value. Since T^ k \uN,0) is negative from the numerical experi- 
ment, liniTv^oo \Tj (ujy, 0)| can be expressed as 

lim 1^(^,0)1 = -^/3cos6<+ r dvsx{v)ln{\ +rj[ k \v)) 

iV->oo 2 J-oo 



oo 
oo 



+ 2 * J f n0 p j l lra l (r)s l{ r) +2\n.thjCi '- (0.1 D 



Finally we obtain the correlation length ( |2.13| ) as 



i -Slnthltf » - f «*»lMln ( • (8.15) 



& 4 i-oo \l + ?7i 1) (^). 



We can solve (|6.5| )-( |6Tl3| ) numerically as follows. First we solve the BAE (3.7) numerically 



for a finite N and determine the Y-functions and their real zeros. This serves as the first 
approximation of their large N limit m and t^ k \ Second we input them into the rhs of 
(|6.5D-(|6Ts|) and get the new rj- functions in the lhs as an output. Third we substitute the output 



^-functions into ( |6 . 1 Q| ) - ( |B . 1 3| ) . Solving them by Newton's method a new output for the zeros 
can also be constructed. Finally by iterating the second and the third processes in the above 
until adequate convergence is achieved, the 77-functions and their real zeros are determined 
accurately. In this way the present approach enables us to overcome the difficulty of the naive 

(k) 

numerical extrapolation of T± as N — ► 00 mentioned in the introduction. 

For a comparison we depict the functions r)^ (v) [k = 1,2,3) in Fig.|| under the same pa- 
rameters. We also include the graphs of the correlation lengths £ 2 and £3 for po = 2, 3, 4, 5 in 
Fig.|5[ In the low temperature limit these results agree with the earlier ones in [12, 25 j. 



/3^oo ZTT& 



with high accuracy. 
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7 Summary and discussion 



We have revisited the thermodynamics of the spin i XXZ model at roots of unity by the 
QTM method. Functional relations indexed by the TS numbers are found among the fusion 
hierarchy of QTM's (T-system) and their certain ratios (y-system). As a peculiar feature of a 



general root of unity, the Y- functions (0)-(4 : 5) and the y-system (|4.6|)- (|4.9|) are considerably 



involved compared with those in [29]. Nevertheless they have a nice analyticity allowing a 
transformation to integral equations. Our approach simplifies the numerics to examine the 
analyticity drastically in that only the largest eigenvalue sector of the QTM T\ is needed for 
the free energy. We have set up Conjecture 1 on the zeros of QTM's supported by an extensive 
numerical study. The resulting integral equations exactly coincide with the TBA equation in 
JTs[ ] based on the string hypothesis. Another and more significant advantage of the present 
method is to allow us to study correlation lengths on an equal footing with the free energy by 
considering other eigensectors of T\. The additional zeros and poles coming into the ANZC 
strips play a fundamental role in characterizing the relevant excited states. We have considered 
the second and the third largest eigenvalues of T\ , which are related to the spin-spin correlation 
lengths for {crja^) and (crjof}, respectively. The excited state TBA equation is derived and 
numerically solved to evaluate the correlation lengths. The result shows a good agreement with 
the earlier one in the low temperature limit. 

Let us remark a few straightforward generalizations of the present results. (1) the XYZ 
model, (2) higher spin cases and (3) inclusion of external field h. For (1) and (2), the T and 
y-systems remain essentially the same. We have an additional periodicity in the real direction 
in the complex u-plane for (1). This does not complicate actual calculations too much. In (2) 
the driving term will enter the TBA equation in a different manner from fl5.7D-(p\9|). As noted 



in [35], the commensurability between the magnitude of the spin and the anisotropy parameter 



would be of issue. This is also an interesting problem in view of the present approach. For case 



(3) the BAE (3/7) should be modified with an extra exp(2/3/i) factor in the rhs. Consequently, 
the BAE roots for the largest eigenvalue will distribute away from the real axis. This is a 
significant difference from the usual row-to-row case where they remain on the real axis even 
for h 7^ 0. The numerical check of the ANZC property therefore needs more elaboration. The 
T-system ( p~^ ) also needs to be modified into 

Ty a+Va _^(v) = T ya ^ ya _^ l (v)+2(-l) mz " C h(f3hy a )T ya _ 1 „ l (v + iy a ). 



Correspondingly, (4.8) is replaced by 

Y ma ^(v) = K{vf + 2di((3hy a )K(v). 

These modifications are consistent with |TJ] from the string hypothesis. Explicit evaluation of 
the effects of the magnetic field on correlation lengths will be an interesting problem manageable 
within the present scheme. 
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Appendix A Takahashi- Suzuki (TS) numbers 

Given {uj} in the continued fraction expansion ( |4.l| ) we define the sequences of numbers 
{^j}^Oi{Pj}fioAVj}j=-iA z j}j=-iA n j}j>ii{^j}j>i and i w j}j>i as follows. The sequence 
{ m j}j=o ^ s define by 

rrij = v\ + v 2 H h z/j < j < a, (A.l) 

m a+ i = oo. 

The sequence }?io * s defined by 

Pj = Pj-2 - Vj-iPj-i 2 < j < a + 1, 

7T 

PO = Pi = 1,P2 = P0 ~ ^1- 

It can be easily shown that 

Pa+1 = 0, (A.2) 
Pj< Et± ,Pj<^- ifl<J<a + l, (A.3) 
2pj + 2p J+ i < p if 2 < j < a or j = 1, i>i > 3. (A.4) 

The sequences {y^ = _ x and {zj}" = _ 1 are defined by 

Vj = Vj-2 + VjVj-x 1 < j < a, 

V-i = 0, yo = 1, yi =^1,2/2 = 1 + ^i^2, (A.5) 
Zj = Zj-2 + VjZj-i 1 < j <a, 

Z-i = l,z = 0,z 1 = l,z 2 = V2- (A.6) 

Obviously, Zj = yj^\\ l , k - tUk+1 and they are all positive integers except y-\ = zq = 0. By 
induction one can verify 

Vj = zjpo + (-l) J p j+ i - 1 < j < a, (A.7) 

Vet = Z a po, (A.8) 

where the latter is a consequence of the former with j = a and (|A.2j) . In fact GCD(y a , z a ) = 1 
is valid. Now we introduce the Takahashi-Suzuki (TS) numbers {nj}j>i fljj and their slight 
rearrangement {%}j>i by 

nj = y r -x + (j ~ m r )y r m r < j < m r+1 , (A.9) 
rij = y r -i + (j - m r )y r m r <j<m r+1 . (A.10) 

Obviously, rij = rij except n mr = y r while n mr = y r -i- In particular, there is a duplication 
n i = n mi = 1) while the modified sequence rij is strictly increasing with j. In this paper we 
are concerned with the first m a + 1 of them. As the set with multiplicity 

{^fHM5 +1 u{j, a }\{i}. 
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We note that if pq > 2 we always have hj = j for j = 1,2 and 3. In parallel with (A. 5) and 
( |A.6| ) we consider the "z-analogue" {iVj}j>i of {nj}j>\' 



Wj = z r _i + (j — m r )z r — 1 m r < j < m r+ i. 



For example, w\ = w 



m-2 



and w 



mi 

n j ~ 1 
PO 



— 1. It is possible to show 
Wj+S jmi l<j<m a , 



(A.ll) 



(A.12) 



where [x] denotes the largest integer not exceeding x. As a result the sequence {wj} is related 
to the parity Vj in (2.14) of @ by vj = (-1)^ for all 1 < j < m a . Using JOD , ( |A~10|) and 
( |A.ll ) one can show 



{lOjPo ± %+i ± 2/r mod 2p } = {po ± (p r - (j + 1 ± 1 - m r )p r+ i) mod 2p } 
for m r < j < m r+ i. Here the signs ± are independent. 



(A.13) 



It is well known IS, 35, 36 that for the TS numbers nj, the equivalent conditions 



,|-U> k ) ) > 0. 



P0 



+ 



raj — k 



Po 



Po 
rij - 1 



Po 



(A.14) 
(A.15) 



hold for k = 1, 2, . . . , nj — 1. It is interesting to observe the condition ( A.14| ) in the light 
of the associated fusion transfer matrix T n -i(u, v) (|3.1|). From (|3.2|) and (|2.5|) , we see that 
( |A.14 ) ensures that e.y and can be independent of their indices for the constituent fusion 
Boltzmann weights to be real. 



Appendix B ANZC property of Y„ 



(i) 



v 



Let us check the applicability of Lemma 1 in section 5 to the Y-system ( [4.6| )-( [0| ) by admitting 
Conjecture 1. Apart from the exceptional Case 1, 2 and 3 listed there, we are to verify that 
Y- (v), for example, is ANZC in the strip S[x] whenever the combination Yj (v+ix)Yj (v—ix) 
takes place. Case 1 has been argued in section 5 and Case 2 and 3 will be considered in Appendix 
C. 

Conjecture 1 tells that the zeros and poles of the Y-functions (|4.3| )-( fL5l) are located as 



y} 1] (v),l + Yj 1 ' (v) : 9to ~ WjPQ ± n j+1 ±y r l<j< 
K^(v): %v ~ w ma -ip ± Ua-l + Vol- 



'(1) 



1- 



Here the signs ± are independent and we have taken the periodicity under v — > v + 2p^i into 
account. See ( |3.8[ ). From (A.13) these functions are ANZC in the following strips: 



For 



m r < j < m r+i — 2 (0 < r < a — 1), 

4 



Y (1) («), 1 + Y (1) ( V ) : 5[p -Pr + (j - m r W - <J], 



1 + (w) : S[p -p r - Pr+l ~ 6], 

1 + Y^v), 1 + kW(v),kW(v) : S[p -p a -5}. 



(B.l) 
(B.2) 
(B.3) 
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In the above 5 denotes a small real number QS\ ~ 10~ 2 ) caused by the deviations of the actual 
zeros of T^2i(v) from the straight line specified in Conjecture 1. They of course depend on the 
Y-functions but have been denoted by the same symbol for the sake of simplicity. On the other 
hand, Lemma 1 is applicable to the Y-system Q4.6| )-(4.9) if the Y-functions are ANZC in the 
strips: 

y} X \v) : S[ Pr+1 ] (m r < j < m r+1 - 2, < r < a- 1), (B.4) 
Y^v) : S[ Pr + Pr+1 ] (1 < r < a - 1), (B.5) 
1 + Y^iv) : S[ Pr - Pr+1 ] (1 < r < a - 1), (B.6) 
1 + Y<£_ 2 (v) : S\p r+ i] (1 < r < a - 1), (B.7) 
1 + *£?(«) : «S[p r ] (1 < r < a - 1), (B.8) 
K^(v) : (B.9) 
1 + K {1) (v), {1 + Yp(v)}™X 1 other than ( B -6)-0 : 5[0+], (B.10) 

where 5[0 + ] means the vicinity along the real axis which can be arbitrarily thin. If v\ = 2, the 
r = 1 case of (|B.7| ) is void. 

Except for Cases 1, 2 and 3 in section 5, it is straightforward to verify that the strips in 
( |B.4| )-( [B.10| ) are narrower than those in ( |B.1[) -( ^3|) for the corresponding functions. As an 
illustration we prove here that S\po — p r + (j — m r )p r+ i — 5] D S[p r+ \] for the strips in (|B.1| ) 
and (|B.4| ). The rest is a similar exercise. We only have to show the inequality 

Po — Pr + (j — fnr)Pr+i — 5 > Pr+i for m r < j < m r+ i — 2. (B.ll) 

Though this is incorrect for j = 1 (hence r = 0), this case corresponds to Case 1, for which the 
difficulty has been cleared in section 5 by a modification of a Y-function. Now suppose j ^ 1. 
It is enough to check (B.ll) only for j = m r (r > 1). Making use of the properties in (A. 3) one 



has 

Po - p r+1 - p r - 5 > po - p 2 - pi - 5 = v x - 1 - 5. (B.12) 



By noting that v\ > 2 and \5\ -C 1, the last quantity is non-negative, proving ( B.ll| ). 



Appendix C ANZC property of Yj (v) for exceptional cases 

Let us show that Lemma 1 can still be applied to the Y-system in Case 2 and 3 in section 5 
after suitable recombinations of the Y-functions. For a function whose logarithmic derivative 
can be Fourier transformed we use the notation 

m(k) = ^- -\nf(v)e- ikv dv. 
2n J-od dv 

We start with Case 2. Explicitly it reads 

yV>[v - i(A- Po ))Y 2 ^(v + i(4- Po ))Y 2 {1) (v-i(p - 2))Y 2 {l \v + i{ Po - 2)) = 

(1 + Yi l \v + i(p - 3)))(1 + Yi l \v - i(po - 3)))(1 + Y 2 (l \v + i(4 - p ))) 

x (1 + Y 2 {l \v - i(4 - po)))(l + Y^\v + »))(! + Y 3 {1) (v - i)), (C.l) 
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where 

r 2 <1 »(t, + i(4-n ) ))r 2 <1 »(»-i(4-p )) 

and the other functions are given by ( |4.13[) and ( (4.14|) . From the Case 2 conditions on V1—V3, 
we have po = 4 — e with < e < 1. Thus the ANZC argument can not be applied to some 
factors in ( |C . 1| ) . For example the (v)-function in the lhs has zeros or poles along Qv ~ 2. 
They are outside of <S[4 — po] Du t can be within S[po — 2]. Similarly zeros of the 1 + Y^(v) 
in the rhs lie along ~ 1 which is in the strip S\po — 3]. A recipe here is to consider the 
combination 

Y 2 ( % + i)y 2 ( %-i) T 3 ( % + i)T 3 ( % - i) 
W 1 + T (v + ^(v - > ^ 

which is ANZC in S[po — 3]. With this, flC.ip can be rewritten as 
X(v + i( Po -3))X(v-i(po-3)) 

= (1 + y 2 (1) (t; + i(po - 4)))(1 + Y^\v - i(po - 4)))(1 + y 3 (1) (t> + i)) 
x(l+y« (<;-*)). (C.4) 

At this stage, the lemma applies to both (CJ5) and (|C.4| ) giving 



F[X\(k) = 2ch^[F 2 (1) ](fc) -T[l+Y^](k), 

Eliminating ^[X] (A;) from these and doing the inverse Fourier transformation, we get (|5[ 
Next we consider Case 3. 

Y} 1] (v + i(l - p 2 ))Y 1 {1 \v - i(l - P 2))y} 1] (v + i(l + p 2 ))Fi (1) (v - i(l + P2)) 

= (1 + y 2 (1) (^ + i))(i + ^ 2 (1) (« - 0)(i + + *(3 - Po))) 

x(l + FW( w -i(3-po))), (C5) 



where 



1 ,y(D M r 2 (1) (z; + ^o-2))r 2 (1) (t;-i(p -2)) 

+ 2 ( j " r ( V + z(3-po))To(,-i(3-po)) ' ( } 



and the other functions are given by fl4T3D and ( gig ). Now p = 2 + p 2 . has zeros 



(resp. poles) at v = ±(1 + u)i (resp. i> = ±(1 — u)i) G <S[1 + p 2 ], which prevents the direct 
application of Lemma 1 in the last two factors in the lhs of (|C,5| ) for J > (resp. J < 0). 



This can be remedied by introducing Y^\v) = y} 1] («)/(th(f "~i|^ u) )th(f " + i|^ u) )) ±jV/2 as in 
section 5. In the rhs of (|C.5| ), there are also some factors possessing zeros or poles and making 
Lemma 1 inapplicable. However the new combinations 

G 1 (v) = (l + Yi 1) (v))(l + Y} 1 \v-i(po-2))) 

T^\v + i(po - 2)) T 2 ( V - i(p - 2)) jfjv - i(po - 1)) 
T^\v - i(3 -po)) r (u + i(4-p )) T (t;-ipo) 
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G 2 (v) = (i + Y ( 1) (v))(l + Y{ 1 \v + i(p -2))) 

T^jv + i(po - 2)) T^jv - i( Po - 2)) T^\v + i( Po - 1)) 
T 1 (1) (v + i(3 -p )) 2b(t;-i(4-po)) T Q (v + ip ) 

are free of these spurious zeros and poles and ANZC in Ssv S [0,1] and [—1,0], respectively. 



With their aid (CI5) can be written as 

YP(V + i(l - p 2 ))^l (1) (^ - *(1 - P2))?1 (1) (V + *(1 + P2))?1 (1) (V " i(l + P2)) 

= Gi(t7 + i)G 2 («-t). 



Solving these relations as in Case 2, we obtain the solution in Fourier space, 



r( i )ir ^ iiVshfcu 
2ch/c 



2cnp2fccn/c 2chp2K 



which can be transformed back to (5.8) 



Appendix D Free fermion case 

Here we consider the free energy and the correlation lengths for the free fermion case A = 
(po = 2, a = 1) in In this case we have <f>(v+4i) = (-l)%(f)(v) and Q(v+4i) = (-l) m Q(v) 



from ( p^q ) and (3.6). Thus (|3.4[) simplifies to 



0(^+2^) 

Ti(u) = g(v) , (D.l 

where 

g(u, v) = 0(v - i(u + 2))0(v + iu) + {-l) m (j)(v + i(u + 2))<f>(v - iu). 

One can directly show T\(v + i)T\(v — i) = (—l) m g(v + i)g(v — i). This rhs is a known function, 
which is a distinct feature of the free fermion model. We find it convenient to introduce 

f[ k \u,v) = 



{v + i(u + 2))<i>(v - i(u + 2)) 

4>(v + iu) + (_jyn 4>(v-iu) \Q{v + 2i) 



It satisfies 



0(« + i(u + 2)) v 7 <f>(y-i(u + 2))J Q(v) 



T[ k \u,v + i)T[ k \u,v -i) = (X(u,v)* + (-l)^~ m X( Ul v)-^) 2 , (D.2) 



where 

_ + i(u - l))(j){v - i{u - 1)) 
[U,V) <t>(y + i(u + l))<l>(v-i(u+l)y 

First we consider the free energy characterized by the largest eigenvalue Ti(u,v). It lies in 
the sector m = N/2. Since the function tP"'(u,v) is ANZC for v E <S[1], we have 

2ch£;.F[f 1 (1) ](A0 = 2T [Jfi + X~\\{k). (D.3) 
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See Appendix C for the notation T. By the inverse Fourier transformation and the identity 

-ikv „ 

-dk = - ; —, = 2tts\(v), 



oo 



e 

oo 2chk 



2ch^f 



we get 

hiTi(u, v) - 
See ( 5.1CQ . Using the relations 

' J (5 



lim X(un,v) 

TV— >oc 



exp 



s 1 * ln(X5 + X~i) (v). 



lim f?°\u N ,Q)= lim T^uj^O), 

N— >oo N—*oo 



(D.4) 
(D.5) 



we obtain the free energy per site / as 

1 



/ 



=- lim lnT^^jv.O) 



2 

7T/9 7o 



In ^2ch(4^ cos 7?) ) (ir/ 



(D.6) 



in agreement with [pfl ]. 

Next we consider the correlation length £2 for (a^a^) which is related to the second largest 



eigenvalue T x (v). This lies in the sector m = N/2 — 1. From a numerical check, T| 
is ANZC for v G S[l]. Therefore we can calculate it in the same way as T^\u,v) The only 
difference is 



u,v 



lnT^O) = 2[si * lnpa - X~*)](v) 



due to (p.2j) with m = N/2 — 1. Thus we have 



2 fi ( 
lim ln\Ti 2 '(u N ,0)\ = - \ In ( 2sh( 

N^oo IT Jo \ 



\JP\ 



cosrj) drj. 



In (2th(fflcos?7) ) d V . 



(D.7) 
(D.8) 

(D.9) 



Combining this with ( |D.6| ) and ( [2.13 ) we find the £2 f° r { a f a i )'• 

1 2 r 
£2 n Jo 

This result coincides with those in Q, |37]-[39|. 

Finally we consider the correlation length £3 for (a^af) characterized by the third largest 
eigenvalue if (u, v). This lies in the sector m = N/2, which is the same with (it, v). All the 



BAE roots for (u, v) are real solutions of g(v) = 0. The set of the BAE roots for T| 
is the same with the one for r x (1) (u, v) except that the largest magnitude ones ±£ are replaced 
by and 2%. It follows from ([Mil) that T^\u,v) = t[ 1) {u,v)th\{v + C)thf (v - (). In the 
Trotter limit, the real zeros ±£ are the largest magnitude solutions to q(un, v)\n—, 00 = given 
by 

-Jf3" 



(3), 



u,v 



c 



-sh" 1 

7T \ 7T 



Thus we have lim jv— «x> l n 
£3 for (cr-of } as 



ifWo) 
1 

6 



lim 



hi 



-21nth -sir 

.2 



2sh 



'(l^l) 



rf^Uiv.O) + 21nth^, and obtain the 

(D.10) 



7T 
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This agrees with |7[]-g|]. 

The results ( p.9| ) and ( p. 10 ) are also plotted in Fig.[|. 
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Figure 2: Location of zeros of T^2i(u,v) for n = 2, 3, 4, 5, po = = 4.8, u = —0.01, N = 16 
(first), iV = 32 (second) and u = 0.01, N = 16 (third). The zeros are located on an almost 
straight line = ±n mod 2p . The deviation from the line is ~ 10 -2 . 
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Figure 3: Location of zeros of T„_ 1 (u, v) (upper) and T„_ 1 (u, v) (lower) for n = 2,3,4,5, 

(2) 

u = —0.1, po = 5, N = 20. T„_ 1 (u, v) has two real zeros for n < 4, which are absent for n = 5. 
T^li{u,v) has two real zeros for n < 4 and a double zero at v = for n = 5. All the other 
zeros are located on an almost straight line = ±n mod 2po. 
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Figure 5: Ratio of the correlation length and the inverse temperature. (<7j~c7~) (upper) and 
(a|c7?} (lower) for p$ = 3,4, 5 and the free fermion case pg = 2. The known result (6.16)-(6.17) 
in the low temperature limit is also depicted by the symbol *. 
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